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The intensities of the 7147 A Ar I and 4806 A Ar II lines obtained from wall-stabilized arcs at pres- 
sures of 0.2 to 5 atm are presented in an Olsen-Richter diagram. Departures from equilibrium are evi- 
dent at electron densities below 5 X 10 Hi cm -3 . The assumption of equilibrium at higher electron densities 
leads to transition probabilities in fair agreement with lifetime measurements. 
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1. Introduction 

For the measurement of optical transition probabil- 
ties and the Stark widths and shifts of spectral lines in 
plasmas no spectral source has been as fruitful as the 
high current atmospheric pressure stabilized arc. And 
of all the elements studied by this technique none has 
been more thoroughly worked over than argon. It is 
therefore a matter of considerable concern that differ- 
ent plasma diagnostic techniques lead to different tran- 
sition probability values in argon arcs [1,2,3].' The two 
plasma diagnostic techniques which have been used to 
establish absolute scales for argon transition probabil- 
ities in arcs are various modifications of the Fowler- 
Milne method [4-7] and the hydrogen line shape 
method [2, 8, 9. 10|. The latter method leads to system- 
atically larger transition probabilities: by about 30 per- 
cent for Ar I lines and by nearly a factor of 2 for Ar II 
lines. Other independent determinations of these tran- 
sition probabilities such as theoretical calculations and 
lifetime measurements generally tend to support the 
higher values for Ar I and the lower values for Ar II. 
Since these discrepancies may be symptomatic of non- 
equilibrium behavior of the arcs it is important to 
discover whether local thermodynamic equilibrium 
(LTE) is a characteristic of these experiments because 
if not, then the large amounts of atomic data being gen- 
erated by the arc technique are of doubtful value. 

If we apply LTE validity criteria estimates [11, 12] 
to the usual experimental conditions existing in argon 
arcs, we find that for electron densities greater than 
5 X 10 16 cm -3 or so we may expect all neutral and singly 
ionized populations to be within at most 10 percent of 
their LTE values. Unfortunately, this LTE threshold 
estimation is uncertain by perhaps half an order of 



1 Figures in brackets indicate the literature references at the end of this paper. 



magnitude because of the approximate nature of the 
theoretical arguments and the rough estimates of 
collision cross sections. Experimental attempts to 
define the LTE threshold electron density have thus 
far not convincingly removed any of this uncertainty. 
The experimental values vary from 5 X 10 14 (Bouras- 
seau, et al.) [13] to 10 1H cm~ :! (Bober and Tankin) 
[6] with most of the measurements suggesting values 
in the range from 2 X 10 ,r> to 10 HJ cm" 5 [14-17]. A 
recent detailed theoretical treatment of the argon arc 
by Uhlenbusch, Fischer, and Hackmann [18] is in 
general accord with the Bober and Tankin value. 

In this paper we present argon arc data which show 
distinct departures from LTE at electron densities 
below 5 X 10 16 cm~ :i and we discuss briefly the evi- 
dence that LTE prevails at higher electron densities 
and the relationship of this problem to hydrogen line 
broadening theory and to argon lifetime measurements. 

2. Experimental Procedure 

Our data consist of spectral intensity measurements 
of two well-isolated, optically thin argon lines. The 
measurements were performed side-on in wall- 
stabilized arcs. The arc chamber for most of the experi- 
ments was of the watercooled copper disk type 
developed by Maecker [19]. Two channel diameters 
were used: 4.8 mm and 3.2 mm. One experiment was 
performed with a 4.8 mm internal diameter fused 
silica tube arc chamber similar in principle to that 
described by Knopp and Cambel [20] in which the 
arc is observed within a transparent constricting tube 
rather than in the open gap between two copper 
constrictors. 

We show in figure 1 our measurements of the inten- 
sities of the spectral lines 7147 A Ar I and 4806 A Ar 
II. At each of the five indicated pressures from 0.2 to 
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4.98 atm experiments were conducted at a number of 
arc currents. Under each operating condition the 
intensities of both spectral lines were recorded at a 
number of wavelengths encompassing the line and at 
generally 40 to 50 equally spaced points traversing the 
arc diameter. Following an Abel inversion to obtain 
local radial volume emission values the line intensities 
were computed as usual by subtracting the continuum, 
integrating over the line, and correcting for the distant 
line wings. Details of the data recording and reduction 
have been described previously [8]. The final result 
consists of corresponding pairs of intensities for the 
two spectral lines at a number of radial points for each 
arc experiment. These intensity pairs are plotted 
logarithmically in figure 1. 

In figure 2 we show curves representing the relation- 
ship between the populations of the upper states of 
these two lines which would exist in complete thermo- 
dynamic equilibrium at the five pressures of our 
experiments. As a matter of interest we have plotted 
also contours of constant electron density and constant 
temperature. These curves were computed in accord- 
ance with Griem's treatment [11] of ionization potential 
lowering, partition function cut-off and pressure cor- 
rection. Since the ratio of a line intensity to the upper 
state number density is a constant proportional to the 
transition probability and since on a logarithmic 
scale such a constant represents a simple translation 
it should be possible to fit the curves of figure 2 to 
the data of figure 1 if the arc data were obtained under 
LTE conditions [5]. The curves in figure 1 show such 
a fitting applied to the data at the three highest pres- 
sures. The fitting was performed numerically by a 
generalized nonlinear least squares fitting technique 



described by Wolberg [21 J. The atomic transition 
probabilities obtained this way are A(7147 Ar l) = 
5.57 ± 0.02 X10 5 s 1 and A(4806 Ar ii) = 7.86 ±0.10 X 
10 7 s -1 . The borders of figure 2 have been selected 
according to these transition probabilities so as to 
permit superposition of figures 1 and 2 and estimation 
of electron density and temperature in the experi- 
mental data. 

Although the least squares fitting was based only 
on the three highest pressures it is clear that the 0.5 
atm data also agree well with this choice of transition 
probabilities. The 0.2 atm data on the other hand are 
in serious disagreement with these transition prob- 
abilities and, in fact, the scatter of data clearly pre- 
cludes any meaningful fitting at this pressure. 

Our data include no corrections for self-absorption. 
Based upon experimentally measured Stark broaden- 
ing parameters [7, 22] and the continuum intensities 
and A values measured here we estimate that the 
maximum error introduced by self-absorption is less 
than 1 percent for the 7147 A Ar I line in our highest 
current, highest pressure experiment and about 1 
percent for the 4806 A Ar II line in our highest temper- 
ature atmospheric pressure experiment. 

By comparing an arc with a black body in which 
complete LTE exists one can recognize three essential 
differences which may cause departures from equi- 
librium. First, in the arc much of the radiation escapes 
so that LTE will be achieved only if radiative depopula- 
tion rates are negligible compared to rates of colli- 
sional population and depopulation. Second, in the 
arc steep temperature and number density gradients 
exist which may generate diffusional fluxes of such 
magnitude that collisional (and radiative) processes 




5 10 lb 2 5 10' * 2 



5 10 ld 2 5 10^ u 2 



4806 flr I I EMISSION 
FIGURE 1. Argon arc line intensity measurements in photons! (cm A s). 

The curves are the family of equilibrium constant pressure contours of figure 2 least squares fitted to the 1.02, 1.98 and 4.98 atm data. 
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4806 A Aril UPPER LEVEL POPULATION (cm -3 ) 

Figure 2. Computed equilibrium upper level populations for conditions of experiments of figure 1. 

Also shown are constant electron density and constant temperature contours. 



are unable to maintain equilibrium distributions. 
And third, since electrons absorb energy from an 
electric field much more rapidly than ions do the 
kinetic temperature of the electron gas will be higher 
than that of the heavy particle gas which obtains its 
energy then mainly by collisions with the electrons. 
Since collisional processes will be dominated by 
electron collisions it turns out that all of these effects 
should decrease rapidly in importance with increasing 
electron density [11, 12]. The importance of the second 
effect should also depend upon arc channel diameter 
and radial position in the arc. 

The data of figure 1 span a range in electron density 
from 10 H; to 5 X 10 17 cm -3 and we see from the fitting 
of the LTE curves that there is a transition from a 
condition where no reasonable fitting is possible at 
low electron densities to a condition where an entirely 
satisfactory fit can be made at higher electron densi- 
ties. This behavior is in accord with the qualitative 
electron density dependence arguments above and the 
assumption that the arc is in LTE at the higher elec- 
tron densities. 

At each pressure the arc currents for the successive 
experiments were chosen to furnish considerable 
overlap between the data obtained near the axis at one 
current and that obtained from the outer regions of 
the next higher current experiment. At the three higher 
pressures and over much of the range of the 0.5 atm 
data there is no significant systematic radially de- 
pendent discrepancy in these overlap regions. The 
data obtained from experiments with the two different 
channel diameters also overlap considerably^ In the 1 
atm data, for example the overlap in the 4806 A Ar II in- 



tensity is two orders of magnitude with the large channel 
providing data from 10 ,r> to 10 19 cm _;i s _1 and the 
small channel from 10 17 to 10 20 cm~ :i s" 1 . The data 
show no evidence that the channel diameter influences 
the results at electron densities above 5 X 10 1(5 cm -3 . 
At lower electron densities we have only the large 
channel data and cannot make the comparison. Our 
results, then are not inconsistent with the equilibrium 
hypothesis at the higher electron densities. The data 
in figure 1 do not include any measurements beyond 
a maximum radial distance which usually lay between 
1.5 and 2 mm. At that distance from the arc axis the 
4806 A Ar II line intensity was generally about two 
orders of magnitude weaker than on axis and the 
remaining data were dominated by noise and scattered 
light. 

The scatter of data in the lower left hand portion of 
figure 1 is caused by a dependence upon the radial 
point at which the measurement is made. This is 
shown in figure 3 where we have replotted the data 
from two typical experiments in this region. Data from 
the arc axis appear in the upper right of this figure and 
data from the outer regions appear at the lower left. 
It will be seen that the deviations from the LTE 
curves fitted at higher electron densities show a sys- 
tematic radial behavior. At the arc center the data lie 
below or to the right of the curve whereas at large 
radii the data lie above or to the left of the curve. 
This behavior agrees qualitatively with that computed 
by Uhlenbusch [18] for a simple model in which the 
kinetic temperature of the heavy particles is allowed 
to be lower than that of the electrons, and the neutral 
atom ground state population is allowed to be larger 
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4806 ArU EMISSION 



Figure 3. Detail from figure 1 . 

Two sub-atmospheric pressure argon arc experiments in 0.48 cm diameter channel: 0.5 
atm, 47.2 amp and 0.2 atm, 62.2 amp. 



than the LTE value obtained from the electron tem- 
perature and from the free electron or excited state 
number densities. 

Using an extension of the Uhlenbusch model in 
which we also include ionic ground state overpopula- 
tions, we can interpret these data in terms of small 
departures from LTE. For this purpose we modify the 
usual equilibrium equations for arc diagnostics which 
relate pressures and temperatures to the populations 
of the excited states and electrons as follows: 
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In these equations the subscript denotes a neutral 
atom; subscript 1 denotes the first ion; and higher ions 
have been neglected for simplicity. Superscript 
denotes a ground state and other superscripts denote 
excited states. n e , m e , and T e are the electron number 
density, mass, and kinetic temperature. Xois the effec- 
tive ionization energy, i.e., the ionization energy minus 
the lowering caused by the plasma interactions. P is 
the pressure and AP is the small pressure correction 
for plasma interactions, bo and b\ are measures of the 
degree of atomic and ionic ground state overpopulation. 
T g is the heavy particle kinetic temperature. We expect 
that T ( ^ T fh bo ^ 1, and b\ ^ 1 with the equality sign 
characterizing the equilibrium condition. 

In order to apply these equations to the data in figure 
3 we introduce the approximation 



and equations for the line intensities 

I$ = A\pJ Q photons (cm ;{ sec) -1 
I{ = A{ n\ photons (cm 3 sec) -1 



(6) 

(7) 
(8) 



= P'-(n l + n e )k(Te-T g )+nl , (bo-l)kT g (5) 



and for the intensity of the continuum which we meas- 
ure near 4806 A: 

/ cont = 1.032 X 10- n w| X-Tr 172 € photons (cm :i s)"" 1 (9) 



where X is the wavelength in angstroms and £ is a 
slowly varying function of T e and X, whose value is of 
the order of magnitude of 1. Now we can use equations 
1 to 9 and the continuum and line intensity measure- 
ments to calculate the ionic ground state overpopula- 
tion factor b\ directly. From our measurements we can- 
not separate the two contributions to P — P' (eq 5) but 
since the two effects operate in opposite directions we 
can obtain a lower bound to bo by assuming T g = T e . Ac- 
tually the assumption T Q =T e probably doesn't introduce 
a serious error near the axis in these experiments. 
Kruger [23] has carried out a calculation similar to 
Uhlenbusch's for a 35 amp atmospheric pressure 
argon arc and finds T e — T g = 200 K at the arc axis 
where the electron density is about 10 16 . Uhlenbusch 
[18] finds a similar value for 200 amperes and an elec- 
tron density of 1.5 X 10 17 . The evaluation of b\ and bo min 
from these equations requires values for the transition 
probabilities and f (jT p , 4806). If LTE is assumed in our 
high pressure data these values can be calculated from 
these data. For this case the A values have been given 
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FIGURE 4. Continuum emission correction factors at 4806 A derived from 1 .02 , 1.98 and 4.98 atm experiments assuming 

equilibrium. 



The curve is a least squares fitted parabola. 



above and values of £ are given in figure 4. The latter 
are in fair agreement with Schluter's theoretical esti- 
mate [24 1 of £ = 2.2 at 14000 K. The results of the cal- 
culation of 6i and b mxn at the axis of the two experiments 
shown in figure 3 are given in table 1. The conclusion 
that, in the 0.5 atm. arc, for example, the neutral ex- 
cited states are underpopulated by about a factor of 2 
and the excited ionic states are further underpopu- 
lated by another factor of 2 is in qualitative agreement 
with the expected non-LTE behavior. 



Table 1 . Departures from equilibrium at axial points of the two 
experiments plotted in figure 3 



p 


T e 


N e 


&o min 


6i 


{atm) 
0.5 

.2 


K 

12300 
14200 


cmr 3 
3.9xl0 lfi 
2.9 X10 16 


1.9 

6.6 


1.9 

4.2 



3. Summary and Conclusions 

With respect to the question of the existence of 
LTE, the presentation of arc spectroscopic data in 
figure 1 can be interpreted to a large extent inde- 
pendently of any assumptions of specific transition 
probability values. Assuming only that the mechanisms 
which lead to departures from LTE become less 
effective as the electron density increases we can 



conclude that there are substantial departures from 
LTE in our low pressure argon arcs. Our high pressure 
arcs are either substantially in LTE or depart from it 
in a manner which consistently mimics LTE over a 
surprisingly wide range of conditions. Although one 
can never expect to obtain complete confirmation of 
thermodynamic equilibrium we feel that the present 
study strongly suggests that arcs in argon at electron 
densities above 5 X 10 HJ cm ~ 3 are in LTE to within 
the few percent scatter of our data. This conclusion is 
based primarily on the apparent disappearance at 
this electron density of the obvious non-LTE effects 
of lower electron densities and the reasonable agree- 
ment of lifetime measurements with the arc-measured 
transition probabilities which result from the assump- 
tion of LTE. The Ar II transition probability measured 
here [A (4806 Ar il) = 7.86 X 10 7 ] is in excellent agree- 
ment with the lifetime-based value (7.9 X 10 7 ) of 
Wiese's critical compilation [3]. The Ar I transition 
probability [A(7147 Ar l)=5.57Xl0 5 s" 1 ] required 
to fit our intensity measurements under the LTE 
assumption is smaller than that indicated by any 
other technique. A comparison of lifetimes for the 
upper level of this line as measured directly [25, 26] 
and as calculated from theory [27-29] and from this 
and other emission measurements [30, 31] is shown 
in table 2. When allowance is made for the possibility 
of systematic errors in the lifetime measurements 
and in the arc measurements — both in the absolute 



75 



intensity measurements of the present experiments 
and in the branching ratio determination — it is doubtful 
that the differences can be considered significant. 



Table 2. An 2p 2 lifetimes 

(nanoseconds) 



Klose [25] 

Verolainen and Osherovich [26]... 

Coates and Gaydon [30] 

Nodwell, Meyer and Jacobson [31] 
Garstang and Van Blerkom [27]... 

Johnston [28] 

Meiners [29] 

Present experiment assuming LTE 



31±2 

28 ±1.4 
*26 

26 

29 

30 

26 
*34 



*Branching ratio taken from arc measurements of ref. [32] 

Many argon transition probabilities have been 
extracted from experiments conducted with a trace of 
added hydrogen which permits the plasma diagnostics 
to be based on electron densities derived from Hp line 
shape measurements. The resulting transition prob- 
abilities, especially for Ar II lines [8, 33], are discordant 
with, for example, the lifetime based values of Wiese's 
tabulation [3J. If the arcs are in LTE then such experi- 
ments indicate that the correct electron density is 
about 10 percent greater than that obtained from the 
Hp line shape measurement and the theoretical 
profiles of Kepple and Griem [34]. This discrepancy 
is roughly equal to the estimated uncertainty of the 
theoretical calculations [34, 35J. On the other hand 
the same kind of contrived departures from LTE 
which would create the illusion of LTE in the data 
of figure 1 but with the wrong transition probabilities 
would also lead to an apparent constant fractional 
error in electron density measurements. Again this 
possibility cannot yet be denied. However, the fact 
that similar apparent errors in ///^-derived electron 
density measurements have been observed in many 
arc gases (hydrogen [36 1, argon [37, 38J, alcohol [39]) 
tend to support the assumption of arc LTE. 
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